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We study a Majorana zero-energy state bound to a hedgehog-like point defect in a topological 
superconductor described by a Bogoliubov-de Gennes (BdG)-Dirac type effective Hamiltonian. We 
first give an explicit wave function of a Majorana state by solving the BdG equation directly, from 
^^ ' which an analytical index can be obtained. Next, by calculating the corresponding topological 

CN ^ index, we show a precise equivalence between both indices to confirm the index theorem. Finally, 

O , 1 we apply this observation to reexamine the role of another topological invariant, i.e., the Chern 

^ ' number associated with the Berry curvature proposed in the study of protected zero modes along 

^^ I the lines of topological classification of insulators and superconductors. We show that the Chern 

number is equivalent to the topological index, implying that it indeed reflects the number of zero- 
energy states. Our theoretical model belongs to the BDI class from the viewpoint of symmetry, 
whereas the spatial dimension d of the system is left arbitrary throughout the paper. 
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Topological classification of matter dates back to reinterpretation of the Kubo formula for the Hall conductance in 
terms of the first Chern number associated with the Berry connection and curvature. [l|, |2| The integer quantum Hall 

^ I effect (IQHE) was thus given a status as a prototype of topologically distinguishable states of matter. An extensive 
use of this idea; reclassifying many-body states by means of topological invariants, has been, on the other hand, 

2 ' started much recently. In the topological classification, matters are classified into a periodic table of ten universality 
classes, including Z and Z2 type topological insulators and superconductors. [S-Q The range of its applicability has 
been further extended beyond such "bulk classifications" , which involves only the Berry connection in momentum 
(fc-) space, to include various kinds of topological defects and associated protected zero-modes. 0, Q 

Q ' With the periodic table in hands, we still need to establish how to relate this topological classification to reality. 

^ [ In the case of IQHE we did not have to ask this question, since a physical quantity, i.e., the Hall conductance was 
directly related to a topological invariant. The same type of reasoning is possible for the Z2 topological insulator, 
-.^ I since it can be regarded as the time reversal counterpart (Z2-version) of the IQHE; as quantum spin Hall insulator. 



J9l4l4j As for the protected zero- modes, similar topological invariant, i.e., the Chern number has been proposed, [7|, 

^\ ' which involves the Berry connection in both k- and x- (real) space. Still, a sufficient reasoning why such topological 

CN| , invariants guarantee the existence of zero modes is lacking. As the Kubo formula for the Hall conductance is to the 

^^ IQHE, so which physical quantity is to the protected zero-modes? We are still in search for its precise analogue. 

l*) [ Specific two models, Dirac fermions with a monopole in 3D [I3 as well as those with a vortex in 2D[la| have been 

[ — ' known as prototypes of fermions with protected zero-energy states. For these models, zero-energy wave functions have 

^^ I been exphcitly obtained. fla [la| and index theorem has been successfully established. jl7Hig | Here, the index theorem 

[20| claims the equivalence between the analytical index concerning the zero-energy states and the topological invariant 

associated with the configuration of the order parameter. The Chern number with respect to the Berry curvature has 

also been calculated for these models. 0, Q Remarkably, the Chern number and the topological index reached the 

same expression, although the starting points are quite different. 

. , For a closer inspection on the proposed Chern number and its precise agreement with the topological index, we 

j^ ' study in this paper a Majorana zero-energy state bound to a topological defect in a BdG-Dirac type superconductor 

in d-spatial dimensions. We first give an explicit wave function of the zero-energy state by solving the BdG equation 

directly. Next, by calculating the corresponding topological index, we confirm the index theorem. Finally, we apply 

this observation to reexamine the role of the Chern number associated with the Berry curvature. We show that the 

Chern number is equivalent to the topological index in arbitrary dimensions, implying that it indeed rcfiects the 

number of protected zero-energy states. 

The structure of the paper is as follows: In Sec. HTJ we construct normalizablc zero-energy wave functions of the 
model, and identify its analytic indices. The equivalence of analytic and topological indices is explicitly established, 
i.e., the index theorem is concretized in Sec. IIIII In Sec. IIVI we reexamine the structure of topological arguments 
based on the Berry connection and curvature, and contrast it with those of earlier sections. We establish an explicit 
one-to-one correspondence between the index of the Hamiltonian and the topological Chern number. In Sec. |Vl 
summary and discussions arc given. 



II. ZERO-ENERGY STATE OF THE BDG-DIRAC HAMILTONIAN 

Let us consider a superconducting condensate described by the following BdG-Dirac type Hamiltonian, 

n = -ir^d, + r''+''cba, (i) 

where j, a = 1, ■ • • ,d. In order to make our parallelism illustrative and convincing The spatial dimension d is kept 
arbitrary throughout the paper. The F-matrices obey {F'^.r''} ~ 26'^'^ for /i, u = !,■■■ ,2d. We also define for 
later use a generalized chiral matrix by F^''+-'^ = (— i)''F^ • • ■ 7^"^, which anti-commutes with F^. The field i^a denotes 
the a-th component of a given multicomponent order parameter. In Eq. ([l} two types of F-matrices, associated 
with the kinetic and the order parameter part of the Hamiltonian, have been introduced in the Nambu notation 
for a system of Dirac fermions in an arbitrary spatial dimension d. When d = 2 the order parameter cf> has two 
components, (p — i4'iy 4'2)y each corresponding to the real and the imaginary part of a gap function. When d = 3 the 
</>-field acquires a third component 03, which represents a mass term encoding a band inversion. Such a 3-component 
order parameter has been proposed for describing proximity effects in a topological insulator/superconductor hybrid 
system. [3, UM The d-component order parameter field 4> in Eq. (IT|) is a natural generalization of this to the case 
of an arbitrary spatial dimension d. In what follows we assume that the chemical potential is set, unless otherwise 
mentioned, to zero; the model has chiral symmetry (see below for details). 

Let us consider a point defect of <p located at the origin, i.e., we assume that (p takes the following hedgehog 
configuration, 

cbjix)=Air)x,, (2) 

where r is the radial coordinate and Xj = Xj/r. We assume A(cx)) = const, and impose A(0) = to ensure the 
regularity of cj) at the origin. Physically, such a point defect may be associated with the presence of a vortex or a 
monopole but here we assume that the hedgehog configuration ^ is given as a background field, and do not take 
account of the effects of associated electromagnetic field. 

The same (mathematical') model, Eq. ([l}, in cases of d = 2 and d ~ 'i, has been extensively investigated, but in a 
different physical context. [iSl - RSj In contrast to these existing studies, the system we have in mind is a topological 
insulator under the proximity effect of a superconductor. [7|, |2l| Let us also mention that a variant of the same 
model with a specific choice of A(r) = r has been recently investigated in Ref. |22|. Compared with this model, our 
model assumes a realistic behavior of the gap function, A(oo) =const., which makes the problem nontrivial even for 
zero-energy states. 

In the classification of topological insulators and superconductors, our setup, i.e., Eq. ([ij with ([2]), belongs to the 
BDI symmetry class, [8| which has both particle-hole and chiral symmetries. Particle-hole symmetry is described by 
{"HjC} = 0, where C is an antiunitary operator with C^ = 1. With an appropriate choice of the F-matrices, it can 
be made C = K, where K is the complex conjugation operator, implying that the Hamiltonian is antisymmetric. 
[23| The chiral symmetry of the model is described by {H,F^''+^} = 0. It may seem awkward to some readers, since 
in conventional relativistic massless Dirac systems, chiral symmetry exists only in even spatial dimensions. In the 
present case, however, since chiral symmetry is extended to include the internal space of the order parameter field; 
dimension is doubled, chiral symmetry is allowed in arbitrary dimensions d. 

We derive normalizable zero-energy wave functions describing electronic states bound to the defect. We deal with 
the case of d: odd and even separately, due to the representation of the F-matrices. 

A. Zero-energy state: case of d = 2n -I- 1 

To solve the eigenvalue equation for H it is convenient to work with a concrete representation of the F-matrices. In 
odd spatial dimensions, d = 2?! + 1, the F-matrices can be chosen as 

F^' = 7^' ® 1 cr\ F''+^' = 1 7^' (g) a^, (3) 

where j = 1, • • • , d, and 

F^'^+l = l(g)l®cr3. (4) 

Here, 7's are conventional 7-matrices given in Appendix|Xl This choice of the F-matrices have the following advantage; 
the Hamiltonian matrix ([1]) decomposes into independent blocks of two chiral components, 

^--f .^W 1. .. ^'^''^ + '^^"f^). (5) 

y 7^ ® lOj - 1 (g) 7-'(/)j J ^ ' 



Correspondingly, let ^ — (\I'_|_, ^_)-^ be a chirally-dccomposcd zero-energy wave function; when 4'„ = (^+ 
the chirality of the wave function ^ is +1 (—1). The eigenvalue equation, W^ = 0, can be written as 

(7-'' ® \dj ±\® 7^0j) *± = 0. 



3 

0), 

(6) 



To solve the above equation, it is convenient to use the following notation: [l5| the action of the matrix a (E> 6 to the 
wave function * is defined as (a ® &*)q,a = o,afib\p^ fip = o-ap'^ pp{b^)p\ = {o-'i'b^)aX, which enables us to regard \l/ 
as a matrix- valued wave function. By the use of this notation, Eq. ([5]) can be written as 



7^9j*±T*±(7^) 0i =0 



If we note Eq. (|A5p . we have 



ya^.^± =P (-)"*±7J' 



0, 



(7) 



(8) 



where we have introduced '^± = ^±r, using F, a product of 7-matrices defined as in (jASp . and n is an integer 
specifying the spatial dimension d ~ 2n + 1. For the hedgehog configuration ^, we may assume that \l/± depends 
only on r, namely, dj'$±{r) = ij^'(r), where \E''(r) = dr'^±{r), and also that the matrix wave function \l/± is diagonal, 

{'i±)ap{r) = Sap'ilj±{r). It then follows that 



Vi(r)T(-)"A(r)V'±(r)=0, 



and its solutions given by 



ip± = C± exp 



±f 



A{r')dr' 



(9) 



(10) 



where C± is a numerical constant. Obviously, only either of V'+ or ip- is normalizable depending on n and on the 
sign of A(c>o). Thus, the normalizable zero-energy wave function is given by 



for (-)"A(oo) < 0, and 



* = 



* 



*^ 







= c. 



C- 



exp 



exp 



(-)" / A{r')dr' 




(-)" / A{r')dr' 



(11) 



(12) 



for (— )"A(oo) > 0. Eqs. (fTTl [T2|) are indeed chirally decomposed zero-energy wave functions, i.e., the eigenvalue of 
Eq. ([4]) is -|-1 for Eq. (fTTj) . whereas —1 for Eq. ([T2|) . Namely, they signify that for a given A(cxd) > the chirality of 
the zero-energy state is (— )"+^. If A(oo) < 0, on the other hand, it is given by (— )". 



B. Zero-energy state: case of d — 2n 



The procedure of constructing the zero-energy wave function is similar. The difference stems from the representation 
of the F-matrices: In the present case, there is not such a convenient chirally-decomposed representation as has been 
adopted in odd dimensions. So that in even spatial dimensions, we use the following F-matrices 



F^ = 7-'' 



1, 



pi+d ^ ^2«+i 



,y, 



where j = 1, • • ■ , d (= 2n), and 



p2d+l ^ ^2n+l 1^ 2n+l_ 



Here, 7-matrices are also defined in Appendix [Xj The zero-energy wave function obeys 

- 17^' (g) 15j* + 7^"+^ ® 7^' (/)_,* = 0. 
By the use of the similar notation for the matrix- valued wave function in the previous subsection, we have 



i7^aj* + (-)"7^"+'*7 



"-^2ji+1^ J 



0, * = *F. 



(13) 

(14) 
(15) 
(16) 



Recall that F is defined in (IA3[) . Taking Eq. ([2|) into aecount and assuming 

$(r) = iir) + i-i'^'^+^ip{r), (17) 

the l.h.s. of the above equation can be converted into 

-i'^i{i(j' + 272"+ V) + (-)"72"+iA(-0 + i72"+V)7-^' 
= -iji [V-' - (-)"Av3] + 7J72»+i [<^' - (-)"A?/'] , (18) 

where -0' = 9r'0 ^nd (^3' = drf. This leads to the following coupled equations 

^'(r) = (-)"A(r)^(r), 

(^'(r) = (-)"A(r)0(r). (19) 



These equations formally allow two independent solutions 



0±(r) = ±(-)>±(r)=C±cxp 



± [ A{r')dr' 
Jo 



(20) 



The normalizability of the wave functions imposes C± = for A(cxd) ^ 0. These two kinds of wavefunctions are 
reminiscent of those in Eqs. pT|) and ([T^. However, in either case of Eq. (PO)) . the chirality of the wave function is 

^ ^2n+lpt5^2n+l 

= (-)"*, (21) 

where we have used the facts that $72"+! = 72"+!$ for $ in Eq. ^ and r1'72"+i = („)n^2n+ipt ^i^^^i^ follows 
from Eq. (|A5|) . We thus find that the chirality of the zero-mode in d = 2n dimensions is given by (— )" independently 
of the sign of A((X)). This is in contrast with the odd dimensional case described in the previous subsection: In even 
dimensions, x — > —x does not cause space reflection and two configurations and —(f> belong to the same topological 
sector. To obtain the zero-energy state with opposite chirality we may consider the following configuration; 

M^) = ^(r) X S^%^ (j = l)^'"'''^^ (22) 

This is a mirror reflection of Eq. ([2]) with respect to the xi-axis. It is easy to see that zero-energy wave function can 
be written as 

* = 7i(V'(r)+Z72"+V(r)), (23) 

where -0 and (p are given by Eq. ((20)) . Due to the presence of 7^ the chirality of the wave function becomes —(—)". 

C. Index of the Hamiltonian: the analytic index 

As in the case of d = 2 and 3, we deflne the index of the Hamiltonian with respect to its zero-energy state. Let N± 
be the number of zero-energy states of H whose chirality, i.e., the eigenvalue of r^'^+i, is ±1, respectively. Then the 
index of Ji is dcflned by 

md'H = N+-N^. (24) 

Let us assume that the Hamiltonian Eq. ([T]) allows only one zero-energy bound state, which has been obtained so far 
in Sec. Ill Al and IHBI Then, provided that A(cx)) > 0, the index ([24]) reads 

■ A-u / (-)" (rf = 2n) ,_. 



III. INDEX THEOREM 

In the last section we have exphcitly constructed chirally decomposed zero-energy states in order to introduce and 
determine the analytical index of the Hamiltonian. The obtained index is given in Eq. (|25p as a function of the spatial 
dimension d. These worked out thanks to specific representations of the F-matrices, e.g., Eqs. ([3]) and @, or Eqs. 
(|13p and (|14p . Here, we consider a priori another index of the Hamiltonian, called a topological index, the evaluation 
of which involves no specific representation of the F-matrices. The equivalence of the two indices is then guaranteed 
by the index theorem: |17H19| 

indT^ = lim TrF^^'+i 



= c+ hm i / .P{x)dS,. (26) 

The first (second) line corresponds to analytical (topological) representation of the index. In the second line, c is 
defined by 

c= lim TtT^^+^—^ ^, (27) 

A/^oo ■H2+M2' ^ ' 

and the current density J^ [x) is generally expressed as 

= i™ ..r-«r. y^^ - ^J^) s'(. - ,). (28, 

where, we have introduced the Pauli-Villars regulator with mass M to make the current well-defined, which yields 
the term c in Eqs. ([26]) and ([27| . [I8| In the second line of Eq. ([26|) the integration is taken over S''"^, i.e., the 
boundary of R''. The first line of Eq. (pS)) . or the definition of index in Eq. ([M)) refiects analytical properties of the 
Hamiltonian, as the solution of the BdG equation as investigated in the previous section. The second line of Eq. (|26p , 
on the contrary, is associated with topological properties of the order parameter. 

The Pauli-Villars regulator introduced in Eq. (|24p becomes particularly important, when the gauge potential is 
taken into account; the second term of (|28|) gives, indeed, the Chern number associated with the magnetic field. 
Here, on the contrary, since the gauge potential is not considered, the integration of the second term vanishes after 
M — > oo. The current density is indeed well-defined even without the Pauli-Villars regulator. (Below, we will verify 
this explicitly.) In the plane wave basis, it can be calculated as 

Here, we have introduced an abbreviation for the inner product, kx = kjXj. In the following, we evaluate this 
expression, i.e., carrying out the integration over the momentum in Eq. (|29|) in two different ways. 

A. Direct calculation of the topological index 

We first evaluate Eq. ([29|) directly. In such a direct calculation, there is no need to specify the representation of 
the F-matrices. Note the relation, 



— ikx'\j2 ^ik 



y^2^^kx ^ _^g^ ^ .^^.)2 ^ ^2 _ ,pr'*+''9, 0,. (30) 



Then, we have 



j3{x) = I -^tTr"'+'p{ir''+-(ba 



(27r)'^ ^ ^^ kj + (I>1+ m2 - ir^r'i+<^d-i(t)a 

f J^trF2^^-lFn^F'^+''^ y. i^T^T^^^^,<f>J 



In Eq. ((3T|) . the series expansion is quite useful, if we note that dj(j)a ~ 0{r~^) as well as dSj oc r'^~^d'^~^^, implying 
that terms with I > d vanish when integrated over the surface of S''^^ in Eq. (|26p . Note also that r^'^+-'^ appears in 
the above equation only once, and hence the terms with I < d ~ 1 vanish if one takes the trace. It follows that only 
the I = d — 1 term gives a finite contribution to the index, which is indeed well-defined. At this point, one can clearly 
see that the Pauli-Villars regulator for the current density, as expressed in Eq. p8|. was indeed not necessary in the 
present case. We reach 

= i'^ r|^j.p2d+lpjpd+apjipd+aipJ2pd+a2 pid-l pd+Cd-ll 

/d'^fc 1 

By the use of tr T^'^+ir^ir^^ . . . p^sd = (2i)''e^i^2--^2'' for ^i^^Vj = 1, • • • , 2d, and 

d-^fc 1 _ r(a-f) 



{2nY{k-^ + hY {An)iT{a)h^ 



(33) 



we have 



(47r)2r(d)((?!)2 +to2)2 

where e^ = (— )''(''~i)/2^ which is due to the rearrangement of the arguments in the totally antisymmetric tensor, 

£i(rf+a)ii(rf+Qi)---(<i+ad-i)jd-i = g^gjji •••jd-i(d+a)(d+ai)-(d+ai_i) ^ Pg^ ^g assume that (j)l -^ |0p = const, at |a;| -^ oo. 
Then, 

0ad, (35) 



2Tl2{d-l)\J\^\ 

(-) edr(^) , aiaa-.-ad. 



„ a,, ,,, , ,!^ai'^V^a2>-'V^a3 • ■ ••^V^ad, 

2tT-2 (d- 1)\ J\x\-^oo 

where (j>a = (j>a/\(f>\, and hence, 0^ — 1. 

So far our calculation was generic, not depending on a specific form of the order parameter, as in Eq. ^. Eq. psp 
is a general expression of topological index for an arbitrary configuration of the order parameter. For a particular 
case of the hedgehog configuration, Eq. ([2|), the integral in Eq. (|35|) reduces simply to the surface area of S"*"^ with 
a unit radius, which leads to 

ind7^ = (-)'^ed. (36) 

When A(c)o) > 0, this coincides exactly with the analytical index in Eq. ((25|) . 



B. Relating the topological index to a transition function 

In the last subsection we have directly evaluated Eq. ([29)) and obtained a couple of expressions for the topological 
index; a general formula, Eq. ([35)) and a concrete expression for the hedgehog configuration, Eq. ([36)) . These 
expressions will be compared with the Chern number due to the Berry connection in the next section IIV Al Here, 
to see a more direct relationship between these two, we derive an alternative expression of the topological index. 
This will correspond to the Chern number in section IIVBI expressed by the transition function. To this end, let us 
introduce another specific representation of the F-matrices than the one we used for the evaluation of the analytical 
indices. The representation we use here makes chiral symmetry of the Hamiltonian manifest. 

Let j^ be the 7-matrices in 2d—l dimensions given in Appendix!^ and let F'^ be the F-matrices in 2d+l dimensions 
which are constructed by 7'' such that 

(M = l,---,2d-l), 

(37) 



F^ , 


= 7'' cr- 


p2d 


= l«)cr^ 


p2d+l 


= 1(g) 0-^. 



Regardless of d, the Hamiltonian turns out to be 



where 



We also introduce a matrix Q by 



^ = ( Pt ^ ) ' (38) 



V = -t-f^dj + ^"^^a - i<Pd. (39) 

^-rkx^^,kx ^Q_^Y.d, (40) 



where V ■ d = T^ dj . It is explicitly given by 

Q = pfcj + r''+°0,. (41) 

For the r matrices Eq. ((37)) we can write it as 

Q^I^Qt^)' Q^ 7'% + 7''+>a - ^0dl. (42) 

The matrix Q plays the role of the transition function in Sec. IIVBI Indeed, by the use of Q'^Q — Rq = J2j=i ^? + 

X]a=i 'Pal ^^ can define the unitary matrix Q = Q/Rq. 

Taking the limit m — !■ 0, we see that the current (|29p can be written as 



where j = 1, • • • ,d. Note that T ■ d = dQ ■ d, where dj = dk , which is distinguished from dj = dx ■ Then, the 
expansion 

= (Q - zr • a)-2 Y.[iidQ-dQ + 2ik-d + a2)-2 (44) 



yields the current 



d'^k 



d 



Ro ^0 \ R'o 



J'{^) = j (^]^tr T^'+'d.QQ^dQ ■ dQY-' + 0{\x\-'). (45) 

Inserting this into the r.h.s. of Eq. (|26p . we can express the index of T-L as an integral over the 2d— 1 dimensional 
space, the d — \ dimensional sphere at the spatial infinities and the d dimensional momentum space. The expression 
thus obtained, however, is not manifestly topologically invariant as it stands. It is possible to write indH as an 
integral of {2d — l)-form as 

where Q — Q/ Rq is the transition function and d = dkjdj + dxjdj is the exterior derivative with respect to x and 
k. The topological invariance now becomes manifest since the {2d — l)-form in the integral changes only by an exact 
form under any continuous deformation of the transition function. For reference, we will give a derivation of Eq. (j46p 
in Appendix [BJ 



IV. TOPOLOGICAL INDICES VS. TOPOLOGICAL INVARIANTS (BERRY CURVATURE) 

How (topological) indices are related to the topological invariants, i.e., Chern numbers expressed in terms of the 
Berry curvature? In Refs. 0, Q, a Chern number, analogous to the one in Rcfs. [l|, |2], but characterizing the 
protected zero-modes, has been proposed. For bulk topological insulators and superconductors, Chern numbers given 
in terms of a Berry curvature in fc-space, are used as a topological quantum number of the ground state. Here, the 
one proposed in Refs. 0, Q involves the Berry connection in both k- and a;-space. 

Let H be the corresponding Hamiltonian proposed in Refs. [3,01, 

H = Pkj + r''+>a + Ar^'^+i = r^x^, (47) 

where fi = 1, 2, • • • , 2d + 1, and X^ = Xfj_{kj,Xj, X). We assume that this model is defined at |a;| — >■ oo, where 
'Pa = l<^l^ = const, so that it includes d parameters {kj}, d— 1 parameters {xj} (these should be denoted by, say, 9j, 
the coordinates of S"*"^, but for simphcity, we use Xj), and one parameter A, i.e, in total 2d parameters {kj,Xj, A}. 
Then, we see that H^ = k'^ + \(j)\'^ + X^ = R^ , which implies that the eigenvalues of this Hamiltonian are ±R. Let ipj 

{j — 1,2,- ■ ■ , 2"^"^) be a set of orthonormal eigenstates with eigenvalue —R, 

H^jjj = ~R^j. (48) 

Then, the Berry connection one-form and curvature two-form are, respectively, defined by 

A = V^d^/>, (49) 

F^dA + A^, (50) 

where ip = {'ipi,ip2, ■ • • ,i^2<'~i) is a 2'' x 2'^~^ matrix and d = dkjdk + dxjdx + dA^A- It should be noted that the 
projection operator [23l to the negative eigenstates can be written as 

P = il^ip'' (51) 

Then, Eq. (ICTI) leads to P(dP)^^P = ^pF'^^p^ which immediately gives 

tr F"^ = tr P(dP)2f (53) 

The d-th Chern number for H is defined by 

^diH)^^(^] ftvF^. (54) 



d\ V27r, 

This is a topological invariant, to be sure, but identifying it as the index is not more than a speculation. As mentioned 
earlier, Eq. (|54p is expressed in terms of a Berry curvature which involves, in addition to kj, the coordinates Xj 
surrounding the defect as parameters for the Berry connection and curvature. This suggests that the Chern number 
([51)1 characterizes the zero-energy states peculiar to the defect. Though in d = 2 and 3, the topological index and 
the Chern number are shown to be equivalent. [19j its direct proof is still missing. Below we establish an explicit 
one-to-one correspondence between the two quantities, which is valid in arbitrary dimensions. The two quantities are 
indeed shown to be equivalent, but possibly differ by a sign factor depending on the spatial dimension. 

A. Direct calculation of the Chern number 

One-form of the projection operator is 



dP = -iP'^dX^, (55) 



where X^ = X^/ R. Then, we have 



trP-^ = -^.irT^X.iv'^dV"' 



^-(2*)V^^--^-+^-i^X^,dX^, ...dX^,,^,. (56) 



22d+l ^ ' P2d+1 



Now let us divide {^^} into {4>a} and {kj} (we regard A as the {d + l)-th momentum kd+i)- Here, note that the 
first Xfi in the above equation cannot be a momentum, since otherwise it vanishes after the integration over the 
momentum. So that 

= ^j^(-)V^"--"-0a,d,/>a,---d0a,d'^+lfc. (57) 

Therefore, we have 



Cd{H) - ^-^ ^^ ^^— ' <:"i"2-a<i 



d\{2nY2d+^ [d - l)\ J ^"-^ ^"-^ ^"^ J (fc2 + |0|2)d+i 

= - o -^^y .M h^'^'-^'k^d^a.---dk.. (58) 

ZTTS (d — 1)! J 

where (/)a = (/)a/|(/'|- It thus turns out tliat 

indH = (-)'*+i+'^(^-i)/2^rf(i?). (59) 

This estabhshes an exphcit one-to-one correspondence between the topological index and the Chcrn number, which 
is valid for a generic configuration of order parameters. The two quantities are indeed equivalent but their relative 
sign shows a systematic alternation as a function of d. For d = 3 this extra sign factor is typically —1. 

B. The transition function 

Here, we give another derivation of Eq. (j59p . i.e., by the use of the transition function. The transition function, 
which will be denoted as Q, includes all the topological information of the Chcrn number. The calculations presented 
in this subsection correspond to those of Sec. IHIBI 

The F-matrices ([57| lead to the Hamiltonian 

where Q is defined in Eq. P^ . Then, Q'^Q + A^l = QQ^ + A^l = R^l. By the use of some reference states, we can 
fix the gauge of the wave functions such that 

^ = P(j)N-'^, (61) 

where N'^ — (f)^ P(j) is a normalization matrix. Let us introduce two kinds of the reference states 

0+ = (?), '/'- = (V)- (62) 

Then, the normalization matrix is given by 



^^± = \/^l (63) 



and corresponding gauge-fixed wave function is given by 

, ^ f -Q 



v/2i?(i? -X)\ Q^ ) 

The implication of the subscripts ± would be clear; ^± have singularities, respectively, in the south and north poles 
on S^'', associated with the Dirac string. It turns out that the transition function ^ at A = defined by 1!}+ = "0-5 is 
given by 

g = i:lip+ = I - Q' (65) 
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where Q^Q ~ 1. Using tr F'^ — A(jJ2d-i{A,F), where uj2d-i denotes the 2d— 1 th Chern-Simons form, and uj2d-i{A + 
V, F) — uj2d-i{^, F) = W2(i-i(V, 0) + da2d-2 with V = dgg~^ and a2d-2 being a certain {2d — 2)-forni, we reach 



In view of Eq. P5| , we have estabhshed the relation ((5^ again. 

V. SUMMARY AND DISCUSSIONS 

We have investigated a Majorana zero-energy state in a topological superconductor in arbitrary dimensions in order 
to reexamine the relation between the topological index and the Chern number associated with the Berry curvature. 
We have first counted the analytical index by solving the BdG equation directly, and next calculated the topological 
index characterizing the topological configuration of the order parameter. We have confirmed that they indeed 
coincide. Based on these observations, we have finally calculated the Chern number also in arbitrary dimensions. 
The obtained Chern number coincides with the topological index apart from an extra sign factor depending on the 
dimension. In this sense, the Chern number for systems with a topological defect is endowed with the role of the 
index for the protected zero-energy states. 

The model we have studied in this paper belongs to class BDI with particle-hole, chiral, and resultant time-reversal 
symmetries. Among them, chiral symmetry has played a crucial role in the index theorem, since the index is defined 
as eigenvalues of the chirality operator. Note that the BDI class with a point defect (in the notation of Ref. H, 
6 = d — D = 1) IS classified by a Z-type topological number, which gives physically the number of protected zero- 
energy states. A more generic class of similar character; possessing particle-hole symmetry with non-trivial topological 
defects is class D, which on the contrary, does not preserve chiral symmetry. In class D protected zero-energy states 
are characterized by a Z2-type topological number, i.e., the choice is simply, whether there exists one of such, or none. 
This reduction of symmetry [25| from class BDI to D is analogous to the relation between IQHE (class A) and Z2 
(quantum spin Hall) insulators (class All) in two spatial dimensions (d = 2, 2D). 

According to the classification by Too and Kane [^ , there are other classes with various kinds of topological defects 
which allow protected zero modes. It may be quite interesting to explore the relationship between the topological 
invariants associated with the Berry connection and curvature, which may be valid, in a sense, in the adiabatic 
approximation, and index-like quantities directly reflecting the protected zero modes, which are free from any approx- 
imations. 
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Appendix A: 7-niatrices 

In this appendix, we summarize the representation of the 7-matrices useful to obtain the zero-energy wave functions 
in Sec. |lTl We first define the conventional 7-matrices in 2n + 1 dimensions such that 



7^^-! = 


= 1®---(8)10O-^ ®fT^ «)■••«) cr^ 


72^= = 


= l(g)---(8)l®o-^(8)cr^®---(g)CT^ 



(Al) 



fc-i 



where /c = 1, • • • , n, and 



n 

which obey {7^,7''} — 26'^'^. Define also, 

r = ^i^3...^2„+i^ (A3) 



and notice 



Then, we have 



for ^= 1,2,--- ,2n,2n+l. 
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(7'^r = (7'^)* = H"r7'^r-\ (A5) 

Appendix B: Derivation of Eq. (|46p 



In the trace of Eq. (|45|) . Q, (9jQ and 9jQ effectively anti-commutc one another because of the presence of r^'^+^. 
Noting this, we can put J^ (x) as 

■^' ("^^ ^ ''j -Bf ^^^' [^''^'^^ ^^^^ ^ ■ • ■ ^^'^ ^^^^^ ^ ■ ■ ■ ^^'^ ^^ 

= (l^/^t^- [r^''^^(d.S)'^S9,.s---a,.Q], (Bi) 

where we have suppressed the nonleading terms for |a;| -^ cxd and dk = dkjdj is the exterior derivative on the 
momentum space. Using this and the volume form on the (d — 1) dimensional sphere 

'i^' = (d^iyi ^""■■''^^^^ • ■ • d^^-^ ' (B2) 

we can write Eq. ((26|) as 

where we have introduced exterior derivatives by d^; = dxj dj and d = da; + d^ . In the representation Eq. p7p of F 
matrices and for Q given by Eq. (J42l). the integrand can be further simplified as 



-l^tr[r2^+iQ(dQ)2'^-i] = -l^tr[QdQt(dQdQt)'^-i _ QtdQ(dQtdQ)'i-i] 

- (-l)'^2tr(QtdQ)2'^-i+d(---), (B4) 

where d(- • • ) denotes total derivative terms and has no contribution to indT^. Eq. (j46p immediately follows from Eqs. 
(1531) and (lB4l). 



[1] D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. den Nijs, Phys. Rev. Lett. 49 (1982) 405. 

[2] M. Kohmoto, Ann. Phys. 160 (1985) 355. 

[3] A. Schnyder, S. Ryu, A. Furusaki, and A. Ludwig, Phys. Rev. B 78 (2008) 195125; AIP Conf. Proc. 1134 (2009) 10. 

[4] A. Kitaev, Proceedings of the L.D. Landau Memorial Conference "Advances in Theoretical Physics", Chernogolovka, 

Moscow region, Russia, 22-26 June 2008 (unpublished). 

[5] M. Zirnbauer, J. Math. Phys. 37, 4986 (1996). 

[6] A. Altland and M. Zirnbauer, Phys. Rev. B 55, 1142 (1997). 



12 

[7] J. C. Y. Teo and C. L. Kane, Phys. Rev. Lett. 104 (2010) 046401. 

[8] J. C. Y. Teo and C. L. Kane, Phys. Rev. B 82 (2010) 115120. 

[9] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95 (2005) 226801. 
[10] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95 (2005) 146802. 
[11] C. Xu and J. E. Moore, Phys. Rev. B 73 (2006) 045322. 
[12] X.-L. Qi, Y.-S. Wu, and S.-C. Zhang, Phys. Rev. B 74 (2006) 085308. 
[13] B. A. Bernevig and S.-C. Zhang, Phys. Rev. Lett. 96 (2006) 106802. 
[14] R. Roy, Phys. Rev. B 79 (2009) 195321. 
[15] R. Jackiw and C. Rebbi, Phys. Rev. D 13 (1976) 3398. 
[16] R. Jackiw and P. Rossi, Nucl. Phys. B190 (1981) 681. 
[17] C. CalUas, Commun. Math. Phys. 62 (1978) 213. 
[18] E. J. Weinberg, Phys. Rev. D 24 (1981) 2669. 
[19] T. Fukui and T. Fujiwara, J. Phys. Soc. Jpn. 79 (2010) 033701. 

[20] For a review of index theorems, see, e.g., G. E. Volovik, The Universe in a Helium Droplet, Oxford University Press, 2003. 
[21] L. Fu and C. L. Kane, Phys. Rev. Lett. 100 (2008) 096407. 
[22] I. F. Herbut and C.-K. Lu, Phys .Rev. B 83 (2011) 125412. 

[23] C. Chamon, R. Jackiw, Y. Nishida, S.-Y. Pi, and L. Santos, Phys. Rev. B 81 (2010) 224515. 
[24] For a recent development of calculations of the Berry connection and curvature, see, e.g., Y. Hatsugai, New J. Phys. 12 

(2010) 065004. See also Ref. [H. 
[25] S. Ryu, A. Schnyder, A. Furusaki, and A. W. W. Ludwig, New J. Phys. 12 (2010) 065010. 



